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The expansion instability of a toroidal current ring in low-beta magnetized plasma is investigated. 
Qualitative agreement is obtained with experiments on spheromak expansion and with essential 
properties of solar coronal mass ejections, unifying the two apparently disparate classes of fast and 
slow coronal mass ejections. 
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We consider the expansion instability of a toroidal cur- 
rent ring with the goal to describe the rapid expansion 
of such rings or partial rings observed in laboratory and 
astrophysical plasmas The equilibrium of this con- 

figuration was established by Shafranov and is realized in 
the tokamak fusion device It necessarily includes an 
external poloidal magnetic field B cx , since the Lorentz- 
sclf force, also referred to as the hoop force, as well as 
the net pressure gradient force of a bent current channel 
always point radially outward. 

The stability of the Shafranov equilibrium has been 
considered by Bateman |4j, who found that the ring is 
unstable against expansion if the external poloidal field 
decreases sufficiently rapidly in the direction of the ma- 
jor torus radius R. Since the hoop force decreases if the 
ring expands, a perturbation dR > will be unstable if 
the opposing Lorentz force due to B ox decreases faster 
with increasing R than the hoop force. Bateman derived 
n = —Rd\nB cx /dR > 3/2 as condition for the instabil- 
ity, which we will refer to as the torus instability (TI). 
The TI can be regarded as a lateral kink instability dis- 
tributed uniformly over the ring. Different from the he- 
lical kink instability, the TI cannot be stabilized by the 
presence of a toroidal field component inside the torus 
(which occurs, for example, in a force- free equilibrium), 
since the hoop force points outward also in this case, with 
an R dependence similar to the purely poloidal configu- 
ration. 

The TI is suppressed in fusion devices by employing ex- 
ternal poloidal fields with sufficiently small decay indices 
n and by stabilizing image currents in the walls of the 
device. However, it may occur in astrophysical plasmas, 
where the external poloidal field is often strongly inho- 
mogeneous Q , and in some plasma experiments 0, 0, • 
In particular, the observations of erupting prominences 
on the Sun, which often evolve into the cores of coronal 
mass ejections (CMEs) causing major perturbations of 
the space weather @, suggest the topology of a single 
expanding partial current ring, whose footpoints are an- 
chored in the inertial visible solar surface. A threshold 
of n > 2 was estimated for this case |J; otherwise the 
instability has apparently never been reconsidered. Re- 
search on CMEs was instead directed at the possibility 



of a catastrophe due to the nonexistence of equilibrium 
in part of parameter space 0, M, ■ 

In the present Letter, we derive a TI threshold that is 
somewhat more general than the one by Batemean and 
treat the evolution of the instability for the first time. 
We consider two cases: a freely expanding ring relevant 
in the laboratory and for CMEs and an expanding ring 
with constant total current, which captures an important 
effect of the footpoint anchoring on an expanding partial 
ring and can be relevant in the initial stage of CMEs. We 
focus on the essence of the instability and its development 
by including only the hoop force (in the large aspect ratio 
approximation, R 3> b) and the stabilizing Lorentz force 
due to B cx . Gravity, pressure, external toroidal fields, 
and any variation in the direction of the minor radius b 
are neglected to permit a largely analytical description. 
The neglect of pressure effects is justified by the fact that 
the instability is primarily relevant for low-beta plasmas, 
in which the conversion of the stored magnetic energy is 
able to drive a large-scale expansion. 

With these assumptions, the force balance is purely in 
the direction of the major radius and given by 0,0] 
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where p m is the mass density of the ring and I is the 
total ring current. The inductance of the ring is given by 
L = noR Qn(8R/b) - 2 + h/2). The internal inductance 
per unit length of the ring Z; is of order unity if the radial 
profile of the current density is not strongly peaked in 
the center of the torus; in particular for uniform current 
density, l\ = 1/2. The flux enclosed by the ring is '5 = 
+ * cx , with 'I'j = LI. Ideal MHD requires * = const 
during a perturbation R — > R + dR. We now have to 
make an assumption how evolves. Here we follow 
Bateman, who ignored changes in the external field due 
to the perturbation and evaluated the enclosed external 
flux using the prescribed external field B ox (R), 

VP = 9 1 + * cx = LI - 2tt / B cx (r)rdr . (2) 

Jo 

This consistency with the use of B ex (R) in the expression 
for the restoring force in Eq. Q implies inconsistency re- 
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garding the conservation of the enclosed flux. If the latter 
were to be treated consistently, one would have to re- 
quire * ex (i?) = const instead. Numerical simulations of 
the instability, which will be reported elsewhere, support 
the instability criterion derived from Eq. (0. They also 
show that magnetic reconnection sets in at the rear side 
of the expanding ring as the instability develops and lets 
the ring effectively "slide" through the external poloidal 
field [lj| , so that Eq. (JSJ) represents a reasonable aproxi- 
mation also for large expansions. With both assumptions 
for ^/ CX (R) it is easily seen that the total ring current, 
I(R) < ^io/L(R), must decrease as a free torus expands, 
since the logarithmic term in L varies only weakly with R 
(subscripts denote initial values here and henceforth). 

We make the ansatz that B CX (R) = BR~ n in the region 
of interest, R> R$. (At R — » a finite B ex is assumed, 
whose particular value will drop out of the equations be- 
low. We also have to assume n ^ 2 in intermediate 
steps of the calculation but find that the final expres- 
sions [right-hand side of Eq. J3J and following] match 
smoothly as (n — 2) — > ±0.) Using Eq. (0), the ring 
current is expressed through the initial values 
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where c = L/(poR). Inserting this in Eq. QJ and nor- 
malizing, p = R/Rq and r = t/T, where 
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is essentially the "hybrid" Alfven time of the minor radius 
(based on the external equilibrium field B eq = _B ex (i?o) 
and the initial density in the torus), we obtain the equa- 
tion describing the evolution of the major radius 
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We now assume c(R) = const, which is exact if the ex- 
pansion is self-similar and can otherwise be expected 
to introduce relatively little error because c depends 
only logarithmically on R/b(R). An approximately self- 
similar evolution of a freely expanding ring has been 
found in a laboratory experiment Q , and also the obser- 
vations of CMEs indicate some degree of self-similarity 
[H|. With c(R) = co, the condition for instability 
d (d 2 p/dT 2 ) /dp\ =1 > becomes 



n > n cx = 3/2 - l/(4c ) . 



(5) 



Bateman's condition is recovered as Co — > oo, which may 
be regarded as the "very large aspect ratio limit." If 
^ ex (R) = $exO is assumed in Eq. (J2J, then the expan- 
sion is described by d 2 p/dr 2 — (c /c) 2 (c + l/2)/(c + 
l/2)p- 2 [l-p 2 - n (c/c )(c + l/2)/(c+l/2)] instead of 
Eq. @ , and [again with c(R) — cq] the threshold rises to 
n > 2. We note that this assumption (with \I/exO = 0) 
and this threshold correspond to the case of a gravita- 
tionally balanced current ring around a star or massive 
object, which should, therefore, be marginally stable. 

Equation J2J can be integrated twice only for small 
displacements, < e = p — l<Cl, showing that the 
expansion starts nearly exponentially, 



e < 1 



(6) 



with the growth rate 7 = (n — n^) 1 ^ 2 . Here vq is the 
initial velocity of the expansion resulting from a pertur- 
bation. Integrating Eq. J3J once shows that for n > 3/2 
a constant asymptotic velocity is reached 



2(2n-3 + ^)(n 



(2n - 3)(n - 1) 



1/2 



[(v T/R f+2] 1/2 , n>3/2. 



(7) 



For n CI < n < 3/2, the acceleration does not decrease suf- 
ficiently rapidly as p — > 00 so that the asymptotic velocity 
diverges. This discrepancy with the behavior at n > 3/2 
results from the simplifications made; it would disappear 
if the restoring forces due to flux and pressure pileup in 
front of the expanding ring, which dominate at large p, 
would be included. The asymptotic gain of kinetic energy 
is AW = M J™(d 2 R/dt 2 )dR ps M(R a /T) 2 , n > 3/2, 
where M = 27r 2 &Qi?op m o is the mass of the torus. 

For a large aspect ratio, the characteristic velocity 
in these expressions is much larger than the hybrid 
Alfven velocity of the initial configuration, Rq/T w 
(R /b )V A i > Vm- Therefore, v T/R < 1 even in 
the case that the initial perturbation vq approaches 
Va; , as may happen if it is due to a kink instability 
The dimensional asymptotic expansion veloc- 
ity « 2 1 / 2 (i? /&o)VAi f° r n > 3/2 scales as the Alfven 
velocity of the initial configuration. 

Figure ^ shows the acceleration profile a = d 2 p/dr 2 
and the numerical solution of Eq. Q with c(R) = Co, 
along with the analytical approximations Eqs. JBJ and 
Q, for particular values of v T/R and Ro/b and for 
the practically relevant range of n. The acceleration rises 
quickly to a maximum, which increases strongly with 
n > n cr and is reached within p < 2 for all n shown. 
It then decreases quickly with increasing p for n > 2 but 
decreases only slowly for n close to n CT . The resulting 
expansion, p(r) — 1, has an approximately exponential- 
to-linear characteristic for n > 2 but is much closer 
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FIG. 1: Radial acceleration profiles and solutions of Eq. Q 
for the freely expanding torus with c(i?) = const and i?o /&o = 
10, voT/Ro = 0.005, ii = 1/2. The approximate solutions, 
Eqs. @ and \Jj, are included as thin lines for p(r). 



to a constant-acceleration curve over a considerable ra- 
dial range for n close to n cr . A qualitatively similar 
n dependence of the acceleration profile is obtained if 
\I>ex(-R) = *exO is assumed in Eq. (J2J. 

This n dependence of the expansion fits perfectly to the 
typical characteristics of CME rise profiles. Fast CMEs 
reach a speed of ~ 10 3 kms -1 , comparable to the Alfven 
velocity in the inner corona, often within a height range 
of h < -R©/3 above the photosphere and show no sig- 
nificant acceleration further out. These events originate 
from active regions which possess a rapid decay of the 
field concentration at heights comparable to the sunspot 
distance D (D ~ i? Q /10 in bigger active regions); for 
essentially bipolar active regions, n > 3/2 for h > D/2, 
quickly approaching n w 3 at h > D. On the other hand, 
slow CMEs propagate with roughly constant, small accel- 
eration throughout the currently observable height range 
(h < 30 Rq), reaching the gravitational escape speed 
of a few 10 2 kms -1 typically only at heights of several 
i?0. These events originate from erupting prominences 
far from active regions, where the large-scale height de- 



pendence of the field, approximately B oc h~ 3 / 2 [l3|. 
dominates already low in the corona. Interestingly, the 
fastest CMEs, and the strongest flares, originate in so- 
called (5-spot regions, which are quadrupolar, with one 
pair of opposite polarity being closely packed within a 
single sunspot, so that a particularly steep field decrease 
(n > 3) occurs low in the corona within very strong 
fields, which imply high Alfven velocities of up to several 
10 3 kms . Thus, the torus instability not only provides 
a uniform description of the apparently disparate classes 
of fast and slow CMEs [l4| but also explains naturally the 
preferred occurrence of the most powerful solar eruptions 
in 5-spot regions fl5| . 

The magnetic field in erupting prominences and CME 
cores can be modeled as a section of a torus, whose re- 
maining part is submerged and frozen in the dense, high- 
beta photospheric and subphotospheric plasma. Such 
line tying is generally regarded to have a stabilizing influ- 
ence; for example, in case of the helical kink instability 
it raises the threshold twist from 2ir to 2.49-7T ^(|. It has 
an even stronger effect on the TI. If a current-carrying 
loop emerges or is formed in the low corona, the line- 
tying is expected to suppress the instability completely 
until the loop is at least semicircular, since the major ra- 
dius of a rising loop must decrease before that stage |l7| . 
Beyond that point, however, the line tying supports the 
expansion because it tends to keep the current through 
the footpoints of the partial ring to be constant. It is not 
clear at present how much of this current can enter the 
coronal part of the ring, where, due to the low resistiv- 
ity, reconnection cannot easily occur within the ring so 
that the number of field line turns and hence IR tend to 
be constant. While a complete account of the line tying 
requires a more sophisticated treatment, we can describe 
its amplifying effect on the expansion by replacing Eq. @ 
with I(R) = Iq, obtaining the limiting case of maximum 
outward acceleration, given by 



dp 
dr 2 
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(8) 



The critical decay index of the external poloidal field for 
instability, 



n cr = 3/2 - l/(2c + 1) , 



(9) 



is only slightly smaller than the critical index for the 
freely expanding ring. The initial evolution is again given 
by Eq. ©. The strong amplifying effect becomes ap- 
parent in the further evolution. This shows an enlarged 
radial range of acceleration, in better agreement with 
some CME observations, and a higher peak (Fig.[5J). The 
asymptotic acceleration does not vanish, however. Since 
a(p — > oo) is small only for ln(8i?o/^o) ^ 1 or for Z i0 ^> 1, 
which both do not have observational support, it is obvi- 
ous that Eq. JHJ cannot hold throughout the expansion. 
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FIG. 2: Radial acceleration profiles of the torus instability 
with fixed ring current [Eq. JHJ] and parameters as in Fig.0 
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FIG. 3: Development of the inverse aspect ratio for the torus 
instability with fixed ring current [Eq. HUH ] for Ro/b a = 10. 



Another consequence of constant ring current is the 
decrease of the aspect ratio in the course of the instability. 
Requiring = * [Eq. witn I(R) = lo, we find 

b(R) _ 8 

R ~ exp{ Co p-i + -p2-n ) + 2 _L} 

n + 2 , (10) 

(Fig. El). Such overexpansion of the minor radius is a 
characteristic of many CMEs, observed as a cavity in the 
outer part of the rising flux, which gives rise to the so- 
called three-part structure of CMEs |8| . The overexpan- 
sion is so rapid that b — > R for p = R/Rq ~ lC^-lO 2 . At 
this point our simplified description breaks down. We can 
expect that magnetic reconnection with the surrounding 
field or between the loop legs is triggered by the overex- 
pansion. This implies that I(R) — const no longer holds 
and that the acceleration falls off as the reconnection pro- 
ceeds. Comparing the acceleration profiles in Figs. and 
12 it is clear that the association of fast and slow CMEs 
with, respectively, high and only slightly supercritical de- 
cay index n holds for line-tied current rings as well. 

Let us finally consider the expansion of a spheromak- 
like torus in a nearly field- free vacuum chamber , which 
proceeded in the observed range p < 2 with roughly 
constant velocity. We note that Taylor relaxation in 
the torus transformed toroidal into poloidal flux in the 
course of the expansion, influencing the TI in as yet un- 
known ways, and that the scatter in the data (Fig. 19 
in Ref. Q) permits a fit with slightly increasing veloc- 
ity as well. With B ox = and ^(R) = LqIq, we obtain 
d 2 p/d,T 2 = (c+l/2)c~ 2 p~ 2 in place of Eq. (@J, where time 
is now normalized to T' = {it / 'co)(&o /Vm) and Vai is de- 
fined using the field (B) at R = and p m o- This acceler- 
ation decreases so rapidly that, soon after onset, the ex- 
pansion velocity is expected to increase only slowly with 
p, consistent with the observation. The asymptotic veloc- 
ity ((c+ l/2)/c 2 ) 1/2 Rq/T ~ 5-16 kms-\ obtained us- 
ing the observed R/b w 2, B ~ 300 G as a representative 
value of the measured range [Figs. 11, 12(b), and 12(c) 
in Ref. 0), and estimated densities N ~ 10 15 -10 16 cm~ 3 



[18j , is in acceptable agreement with the observed expan- 
sion velocity of w 5 kms . 

We conclude that the TI is a possible mechanism for 
CMEs (in addition to a catastrophe 0, and to the 
helical kink instability Jicj), that the TI governs their 
medium-scale (p < 10 2 ) expansion, providing a unified 
description of fast and slow CMEs and a possible ex- 
planation for their three-part structure, and that the TI 
occurred in experiments on spheromak expansion. 
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